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1. Introduction 

Current observational data support standard ACDM model greatly [1], in which primordial perturbations which 
are assumed to be responsible for Cosmic Microwave Background anisotropics and large-scale structure for- 
mation are generated from quantum fluctuations and stretched to superhorizon scales during inflation (see e.g. 
[2] for a review). Standard single-field slow-roll inflationary models predict almost scale-invariant, Gaussian 
and adiabatic primordial fluctuation. 

Actually, inflation itself is not a single model, but rather a theoretical framework. From the point of view 
of power spectrum, many inflation models are "degenerate". Power spectrum, or strictly speaking, tree-level 
two-point correlation functions of cosmological perturbation do not give us an unique theory of inflation. In 
fact, phenomena beyond linear-order have been extensively investigated over the past several years. 

The most significant progress beyond power spectrum is the investigation of statistical non-gaussianities of 
CMB anisotropics and primordial fluctuations (see e.g. [4] for a nice review). From the field theoretical point 
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of view, non-gaussianity describes the interactions of perturbations, which will cause non-vanishing higher- 
order correlation functions. Such interactions are mandatory in any realistic inflationary models, which come 
from both the non-linear nature of gravity and the self-interactions of inflation field(s). In standard slow- 
roll inflation scenario, however, the non-gaussianities have been proved too small to be detected [3], even 
with PLANCK. Thus, any detection of non-gaussianties would not only rule out the simplest inflation models 
but also give us valuable insight of fundamental physics of inflation. Various models have been investigated 
to generate large non-gaussianties by introducing more compUcated kinetic terms [5, 6, 7, 8] or more fields 
[9, 10, 11, 12, 13, 14, 15, 16, 18, 19, 22, 23]. 

Interactions not only cause non- vanishing higher-order correlation functions which are responsible for non- 
gaussianities of primordial perturbations, they also cause quantum loop-corrections. It is possible that such loop 
corrections may be large in their own right. Especially, we could expect large loop corrections in models which 
have also large non-gaussianties, since both of them describe the interactions among perturbations. Moreover, 
by collecting signatures of both quantum loop corrections and non-gaussianties, we will obtain a more sensitive 
test of the physics of inflation. 

Early estimations of loop effects in cosmological perturbations were pioneered by Mukhanov, Abramo and 
Brandenberger [24, 25], Abramo and Woodard [26], and Unruh [27]. Subsequently, loop effects have been in- 
vestigated by many authors [28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49]. 
Recently, stimulated by works of S. Weinberg [50, 51], loop corrections to primordial fluctuations have been 
re-investigated by several authors [52, 53, 54, 55, 56, 57, 58, 59]. It has been found that cosmological loop 
corrections seem to have some infrared divergences, which scale as ln(A;^), where is some IR comov- 
ing momentum cut-off [52, 53, 54, 55]. Other discussions on cosmological loop corrections and the infrared 
divergences can be found in [63, 64, 65, 66, 67, 68, 69]. The presence of IR divergences implies that loop- 
contributions may potentially give large corrections to tree-level results, although it was argued in [63, 64, 65] 
that these IR divergent corrections do not contribute in quantities that are directly observable. 

Actually, now there are two types of so-called "loop corrections" in the literatures. One type is the quantum 
loops (q-loop) which originate from the interactions among fields and are evaluated using standard quantum 
field theory, while the other type is the classical loops (c-loop) which arise mathematically from some "loop- 
type integrals" due to the non-linear map from Scj) to C, [60, 61, 62]. In this work, we focus on the quantum 
loops. However, it is well-known that the observational variable is the curvature perturbation ( and there 
has non-linear relation between ( and 6 during superhorizon evolution (e.g. in terms of the "(^A^-formalism" 
[73, 74, 75]). Thus, these "loop" effects also have to be treated seriously ,and especially they can be significant, 
if the higher-order derivatives of e-folding number N with respect to the inflaton fields are extremely large. 
While in [62] it was shown that the loop corrections for ( are always suppressed during inflation, if one assume 
that the fields follow a classical background trajectory. 

Up to now, one-loop corrections to the power spectrum of the cosmological perturbations have been inves- 
tigated in single scalar-field loops [52, 53, 54, 55, 50, 51], concerning effects from many light scalar fields [57], 
and from graviton loops [56]. Actually, it has been recognized that it is hard to generate large "local-type" non- 
gaussianities in single-field models, unless we abandon the slow-roll-type conditions. Thus, the investigation 
of multiple field models has particular importance. In multiple field inflationary models, perturbations can be 
decomposed instantaneously into one adiabatic mode and several entropy modes [20] (see [21] for a review). 
There are interactions between adiabatic mode and entropy modes. 

In this work, we investigate the one-loop corrections to adiabatic power spectrum due to the interactions 
between entropy mode and adiabatic mode, at the horizon exiting. We find that due to the non-canonical kinetic 
term in DBI models, the loop corrections are enhanced by slow-varying parameter e and small sound speed 
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Cs, as the enhancement of non-gaussianties in such models. Thus, in general, the loop-corrections in multi- 
DBI models can be large. Moreover, we find that the loop-corrections from adiabatic/entropy cross-interaction 
vertices are IR finite. This is because in the limit of speed sound c,, the leading-order three-point interactions 
are totally "derivatively interacted". While the derivatives (both with respect to time and space) give momentum 
factors which cancel the momentum factors in the denominator from the tree-level correlation functions to make 
the loop-integrals IR safe. 

Quantum cosmological correlation functions are evaluated by using "in-in formalism" (also dubbed as 
"Schwinger-Keldysh formalism" ) [70, 71, 72] (also see e.g. [50] for a nice review), which we also give a brief 
review in Appendix A. In [58] the application of "in-in" formalism on cosmological perturbation, especially on 
loop corrections, has been investigated in details. 

This paper is organized as follows. In section 2, firstly we briefly review the linear perturbation of multiple- 
DBI models, then we describe the third-order perturbation action, especially the three-point adiabatic/entropy 
modes cross-interaction vertices. In section 3, we calculate the one-loop corrections to the adiabatic power 
spectrum from the adiabatic/entropy modes cross-interaction vertices, by using the "in-in" formalism. Finally, 
we make a conclusion and discuss the related issues. 

In this paper, we set = c = SttG = 1. 

2. Basic Setup 

Perturbation theory of general multiple field modes with non-canonical kinetic terms up to third-order action 
has been calculated by several authors [10, 11, 12, 13, 14, 15, 16, 17, 22]. In this work, we focus on multi-field 
DBI models [10, 11, 12, 13, 14, 15, 16]. In this section, we briefly review previous known results and setup our 
conventions. 

2.1 The model and the background 

In this work, we focus on two-field DBI inflationary models, with action of the form 

d^x i\R + cA , (2.1) 



with! 



where 



Cf = -^7=y ( v^l - 2fX + 2/2 (X2 - XijXiJ)^ - V{cj>') , (2.2) 



X'^=--d^(p'd^^(p-\ (2.3) 



andX = GijX^-^. 

The slow-rolling parameter is defined as usual 



H 



(2.4) 



'in general, the full DBI-type action takes the form C = —j (^Vd — Ij —V, with 

D = det {Si J - 2fXij) = l-2fX + 4fxl'X;^ - 8fx\'xjX^^ + IsfxYxjX^X^'^ . 

However, it is shown that higher-order terms proportional to and /■* do not contribute to the calculation of power spectra and 
leading order bispectra. Thus, in this work we neglect these terms for simplicity. One can refer [10, 11, 12, 13, 14, 15, 16] for detailed 
investigations of multi-DBl inflation models. 
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The background "inflaton velocity" a can be related with e as 



a = ^riX = H^^s, (2.5) 
where Cg plays the role of the propagating speed of the perturbation. In DBI inflation model, 

= 1 - 2/X , (2.6) 
where X is evaluated with the background quantities. 
2.2 Brief review of linear perturbation 

In multi-field inflation models, perturbations can be instantaneously decomposed into one adiabatic mode and 
several entropy modes [20]. Restricting to two-field case, we introduce the basis {6^,65} with = jt, 
where a is the background inflaton velocity defined in (2.5), and e\ is orthogonal to e^. Then we introduce the 
decomposition 

= S^' = Q^ei + Qse\. 7=1,2. (2.7) 

Qcr and Qg are adiabatic mode and entropy mode respectively. 

After the adiabatic/entropy modes decomposition, the second-order perturbation action for two-field DBI 
inflation model up to leading-order was derived in [10, 11, 14] in spatially-flat gauge: 

5(2) = I d^d^x ^ { [(Q;)2 - cl{d,Q,f\ + cl {{Q'gf - clmsf\ } , (2.8) 

where if = d/a is the Hubble parameter, r\ is the comoving time, and a prime denotes the derivative with 
respect to t]. In deriving (2.8), we assume that is approximately constant and neglect coupling between 
adiabatic and entropy modes which may give significant effect in some cases, and keep only the leading order 
terms. Note that in general multi-field models, adiabatic mode and entropy modes can propagate with different 
(and arbitrary) speeds of sound, which was first pointed out in [12, 13] (see also [14, 15, 16, 17, 22]). Multi- 
field DBI model is a special case, where adiabatic mode and entropy modes propagate with the same speed of 
sound C5. 

In canonical quantization procedure, we write 



Qaij], k) = Uk{r])ak + ul{r])a^_f^ , 
Qsiv, k) = Vk{ri)ak + vl{r])a'_f^ . 

In de Sitter approximation a{r)) = —l/Hrj, the mode solutions can be easily solved 

Vkiv) = /— 3 (1 + icskv) e-'"''''^ . 



(2.9) 



(2.10) 



Note that this choice of mode functions corresponds to the Bunch-Davis vacuum. Now it is straightforward to 
write down the tree-level two-point correlation functions (see Appendix A for more details), 

{0\Q^{v,ki)Q^{r]',k2)\0) = {27Tf5''{ki + k2)GkAv,v'), 
{0\Qs{v,ki)Qs{v',k2)\0) = {27rf6^{k, + k2)Fk,{v,v') , 



-4- 



where 

Gk{ri,ri') = Uk{v)ul{v) , Fk{v,v) = Vk{v)v*k{ri') , (2.12) 

are sometimes called Wightinan functions. 
Then the power specti-a for Q^r and Qs are 



Pa{k) = \Qa*\^ 



Hi 

Hi 
2k^d ' 



2 (2.13) 



respectively, where the subscript * indicates that the corresponding quantity is evaluated at the sound horizon 
crossing kcg = aH. 

2.3 Adiabatic-entropy modes three-point cross-interactions 

The third-order perturbation action of two-field DBl models was derived in [10, 11] (see also [12, 13, 14]). In 
the leading-order of slow-varying parameters and in the limit of small c^, we have 

The first line in (2.8) is Qcr three-point self-interaction terms, while the second line is Qa-Qs cross-interaction 
terms. Note that there is no Qs self-interaction term. 

In this work, we focus on the cross-interaction terms involving entropy mode and adiabatic mode. The 
interacting Hamiltonian at third-order is given by Hi = — ^(3), where L(3) is the Lagrangian read from (2.14). 
There are three cross-interaction vertices, 

Hi = Hia + Hib + Hi,, (2.15) 
where subscript "i" denotes interacting part of the Hamiltonian, and 

^-(^) = f /n^(27r)=^5'(fci23)QU^>fci)Q'.(^,fc2)Q',(r?,fe3), 

Hibiv) = / n ^ (27r)^<5=^(fci23) (fe2 • ks) Q'.ir,, k,) Q,{ri, k-^) Qs{ri, k^) , (2.16) 

HUV) = ^ / n ^ {2T,f6'{ki2s) (fei ■ fc2) Qa{ri, fei) Qsiv, ^2) Q'siv, ks) , 
where we have changed into momentum space and 



^In the literatures on "in-in formalism", it used to define "two" types of Green's functions: G^(?7, r;') = Wfe(??i)M*(»72) and 
Gkiv^v') = Mfc('7i)w(»?2), keeping r? > r]'. It is easy to see that, mathematically, G^{rj,r}') = (0^(77,77'))* = 0^(77', 77). Thus, 
if we relax the artificial restriction 77 > 77', there is no need to keep track of this odd notation. In this paper, we use Gk{r]j v') for 
simplicity. 
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which plays the role of an "effective coupling constant" and we assume to be approximately constant during 
inflation. In deriving (2.16) we have used de Sitter approximation a = —l/Hrj, and fci23 = fei + ^2 + ^3 for 
simplicity. 

From (2.16) it can be seen that the three types of three-point interaction vertices are all "derivative in- 
teractions", as shown in Fig. 1. This is different from ordinary field theory where fields are locally coupled. 
Moreover, as we will see, the derivatives bring momentum factors, which will cancel the momentum factors in 
the denominator come from the two-point correlation functions, and increase the finiteness in IR. Especially, in 
this work we will show that loop corrections arising from the above three-point vertices are all IR finite. 

/ a _/ a 



fir \ G,. \ 6^ 

Q, a 

Figure 1: Diagrammatic representations of three types of adiabatic-entropy three-point cross- 
interaction vertices. A red dot denotes temporal derivative while a blue dot denotes spatial deriva- 
tive, or in fourier space, a momentum factor. 

There is another issue should be emphasized. In standard slow-roll inflation models, the "effective cou- 
plings" of the interaction vertices are greatly suppressed by slow-roll parameters (this is also why primordial 
non-gaussianies are much small in these models) and thus contributions to the correlation functions from multi- 
vertices diagrams are suppressed and can be neglected at least in the leading-order. However, this is no longer 
the case for Cg <^ I models, where the strength of the interactions are much larger. This can be seen directly 
from (2.17), the effective coupling is enhanced by small Cs- For example, in [23], trispectrum from exchanging 
scalar modes were investigated, it was shown that in ^ 1 models, in general the contribution from two- 
vertices diagrams cannot be neglected in comparison with one-vertices diagrams. The essential idea is that, 
contributions from diagrams with more numbers of vertices are of the same importance in comparison with 
contributions from diagrams with less numbers of vertices, even in the leading-order. This is also one of the 
motivations of this work. 




3. One-loop Entropy Mode Corrections 

Correlation functions in cosmological context are calculated by using "in-in" formalism [70, 71, 72] (see Ap- 
pendix A for a brief review). In this work, we focus on the one-loop contributions arising from the three-point 
vertices, thus we need 



{Qa{v*,f^i)Q<T{v*,k2))i-ioop = -2Re 



V* rvi 

dvi / dr]2 {0\Qa{ri^,ki)Q^{r]^,k2)Hi{rii)Hi{r]2)\0) 

-00+ J —00 + 

rn* rv* 

+ drill dr]2 {0\Hi{r]i)Q^{r]^,ki)Q^{r]^,k2) Hi{r]2)\0) , 

J —oo~ J —00+ 

(3.1) 

where we have written explicitly the lower bound for time-integrals^: — oo^ = — cxd(1 =p ie) (see Appendix A 
for details). Note that all quantities on the right-hand side in (3.1) are "interaction-picture" quantities. In this 
paper we choose Hi involving three-point cross-interactions of Qcr and Qs which is given by (2.15) and (2.16). 



^One can refer a recent work [58], for a detailed discussion on the application of "in-in" formalism on cosmological correlation 
functions, and the issue of the lower bound of time integrals. 



-6- 



In the following, we neglect the subscript "i" for simpUcity. Denote 



V* rvi 



{Qa{ri*,ki)Qa{v*,k2))af3i = -^^^ I dr]i I dr]2 {Q\Qa{v*,ki)Qa{v*,k2) Haim) Hp{r]2)\Q) 

\_J —00+ J —00+ 

/V* rv* 
dm I dr]2 {0\Ha{m) Q<7{m,ki)Qa{v*,k2) Hi3{r]2)\0) , 
-oo~ J —00+ 

(3.2) 

where i^a's are given by (2.15) and (2.16). Thus, the final loop-corrections read, 

2 

((3a(??*,fci)Qa(??*,fe2))i_ioop = XI ^ {Q<'i'n*,ki)Q„{r]^,k2)) ■ (3.3) 

a,P=a,b,c 1=1 

It is useful to note that (• • •)^^ = (• • •)^^. 
3.1 Calculating the Loop Corrections 

Now we investigate the loop corrections. The calculation is straightforward but rather lengthy. In this section 
we show the explicit steps. Readers who are interested only in the final loop contributions can go to the next 
subsection directly. 

3.1.1 Diagonal contributions 

From (3.1), the first integral reads. 



{Qa{V*,ki)Qa{v*,k2))aal = "^Re 



/V* rm 
dm / dm {o\Qaim,ki)Qaim,k2)Ha{m)Ha{m)\o) 
-oo J —oo 

(3.4) 

The integral in the squared-bracket in (3.4) is 



/n* rm Id d 

dm / dr]2 — -T-GkAr]*.m) -T-GkA^*,m 



) 



f d^p ^ , , d^ ^ 

J (2^^'^^'^''^^H^'^^^-'^i^'^''^^- 



(3.5) 



Actually, this expression can be read from the "Feynman-type" diagramatic representation of the contribution, 
as shown in Fig.2. The main differences from usual field theory is that, here the perturbations are mostly 
"derivatively-interacted", rather than in usual field theory the interactions are described by local products of 
fields. Moreover, in cosmological context, we are interested in the expectation values of some arbitrary time 
rather than scattering amplitudes. This is why the so-called "in-in" formalism has to be taken. The only thing 
we should note is the order among three "time" — ry*, m and r]2 between (0| and |0). More precisely, in 
calculating the first integral of (3.2), the order is (from left to right) 'q^-m-V2, while in the second integral in 
(3.2), the order is m-V*-il2- Keeping these things in mind, we are able to write the above exphcit expressions 
from Feynman-type diagrams directly. 
Now, from (2.10) and (2.12), we have 

:^GkAv*,m):^GkAv*,m) = ^ vm {i + icshv*f e"^^'''^'''+^'-^'^*\ 



dmdr]2 d?7id?72 
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p 



p-ki 



Figure 2: Diagrammatic representation of the (aa) loop contributions. A "red dot" 
denotes derivative with respect to comoving time. It is important to note that the left 
and right end-points are labeled with 77*, while the left interaction vertex is labeled 
with 771 and the right vertex is labeled with 772 ■ 

Thus, (3.5) can be written as 

(27rf 5\k^ + k2)gl^il + icshv*f e-^-^'^i^* J -0^ p \p - k,\ 

-00 J —00 



(3.7) 



The lower bound of both time integrals are —00, where the integrants highly oscillate. However, we have to 
remind ourselves that the lower bound should be chosen carefully in order to choose the true vacuum (see 
Appendix A for a detailed discussion on this problem and the issue of "choosing vacuum"). Actually, here the 
lower bound of rji integral and rj2 integral is the same, —00+ = —00(1 — ie). Thus, the time integrals can be 
evaluated without any ambiguity. For example, the integral w.r.t. 7]2 is evaluated as 

,1 , + ; gg'-a , (3.8) 

— oo(l— ie) 

where we write K = ki+p+\p — k\\ior: simplicity. Thus, after performing the time integrals and taking the 
real part, we finally get 



(Q.(r?*,fci)g.(r?.,fc2)U = -(2vr)3 53(fc^ + fc,)^2^ j ^P\P_p[ (a,K^ + a2KK + a,k^ 



c^H^ f d^p p\p — ki\ 

'a 

where we have defined K = ki — p — \p — ki\, and 



(3.9) 



ai = 2 (2x1 - 5a;2 - 5) , 



a2 = -{2xt + bxl + b) , (3.10) 
as = - {2xt + + 1) , 

with X* = —Cgkim for short, note that K + K = 2ki. 

The next step is to evaluate the momentum loop-integral: 

I{k) = j ^'^^3^'' {aiK^ + a2Kk + 03^') , (3.11) 

remember that K = ki+p + \p — ki\ and K = ki—p—\p — ki\. Note that the integral (3.11) is IR finite, thus 
we do not need to introduce infrared cutoff here. It will be shown in this work that loop corrections from 
the three-point vertices from (2. 16) are all IR finite, due to the "derivatively-coupling" of the perturbations. The 
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integrant in (3.11) has the structure f{ki,p, \p — fei|), in Appendix B, we describe the method to evaluate such 
type of integrals. In particular, here the function f{ki,p, |p — fci|) is 

f{ki,p, \p-ki\) = (aiK^ + a2KK + as^') • 

After a straightforward but tedious calculation, we get 

(Qa(r?*,fci)Qa(r?*,A:2)>a„i = {2TTf6''{ki + k2)Pa{ki)^^{cilnk + ai) , (3.12) 

TT eCg 

where 



" 1920 ' 

and ai denotes the finite part of the momentum integral together with the scheme-dependent UV renormaliza- 
tion constant, which can be found in Appendix C . 

Now let us move to the calculation of (aa2) contribution. The expression can also be read from fig.2: 

{QAkl,V*)Qa{k2,V*))aa2 

/V* rv* I r ^ ^ 

-J"' J- J''' ^ y d;^^'^^^^^'^*) ^2^'^^"*'"'^ (3.14) 

d^ . . d2 



drjidr]2 dr)idr)2 

It is useful to note that, the only difference of integrant from (aal)-contribution is that, here in fig.2, the order 
of time is rji — rj* — r)2. That is, what we need to do is to simply change the first Green's function in (aal)-case 
as following 

-7—GkAv*,m) =^ -7—Gki{vi,v*), 

drii dr?i 
to get the above result. Following the same strategy as above, we finally get 

(Qa(fei,7?*)ga(fe2,??*))„„2 = {2TrfS\ki + k2)Pa{ki)^^{c2lnki+a2) , (3.15) 

with 



C2 = ^xt{l + xl) , (3.16) 

and a2 can be found in Appendix C. 

Collecting (3.12) and (3.15) together, we get the contribution from (aa)-term 

{Qaiki,ri,)Qa{k2,V*))aa = {27r f {ki + k2) Pa (kl) (Caalnki+aaa) , (3.17) 

TT eCg 

with Pa{ki) is the (tree-level) power spectrum of adiabatic mode Q^, and 

Caa = Ci + C2 = ^ (-2 - 2x1 " Wxj + 3xt) . (3.18) 
The contributions from (bb)-term can be read from Fig.3. 

{Qa{kl,V*)Qa{k2,V*))bbl 

= -{27rfdHk, + k2) Sg^4 X Re£" j^'' dm ^ j ^ [p • (p - k^)f (3.19) 

X -^Gk^{r]^,r]i) -^GkA'n*,m) Fp{ni,m) F\p-ki\{m,m) ■ 

QTji ar]2 
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p 



' ' 1 11 m\» • 



p-ki 

Figure 3: Diagrammatic representation of the loop corrections from (bb)-term. As 
before, a "red dot" denotes derivative with respect to comoving time associated with 
the vertex, a "blue dot" denotes a momentum factor associated with the line mo- 
mentum. We may freely label the momentum flows into a vertex a "+" sign and 
momentum flows out of a vertex a "— " sign. The momentum factors of the same 
vertex are "scalar-producted", e.g., the left vertex gives a factor p ■ {p — ki). 



And 



{QAki,V*)Qa{k2,7]*))bb2 = i'^^f^^ik, + k2)Ag^ct r dm r dm-^ [ -^\p-{p-ki)f 



Aril d7]2 

(3.20) 



After lengthy but straightforward calculation, we get 



jj2 

{QAki,m)Q^{k2,v*))bb = {27rf 6^ {ki + k2) Paih)^^ {cbblnki + abb) (3.21) 



with 

-273 + 27x1 + 90xf - 8x5 



Cbb 



3840 



(3.22) 



agian abb is the finite part of the momentum integral and can be found in Appendix C. 
Contributions from (cc)-term can be read from Fig.4. 



ki 4 ^ k2 ki 4 ^ k2 



p-k, p-k] 



Figure 4: Diagrammatic representation of the loop corrections from (cc)-contribution. As before, 
a red dot denotes derivative with respect to comoving time associated with the vertex, a blue dot 
denotes a momentum factor associated with the line momentum. Note that there are two types of 
diagrams from (cc)-contributions. 
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((3^(fcl,7?*)(3^(fc2,^*))ccl 

/V* rvi I r d^j) 
dvi / dv2 / TTT^ Gk, {v*,m)Gk,{r],,r]2) 
-oo J-oo mV2 J (2vr)'^ 



d2 (3.23) 



r 2 

-{p-ki) [fci-(p-fei)] ^^p(^i>%) ^^p-fci|(f?i,%)| 



And 



{Qa{kl,'n*)Qa{k2,V*))cc2 

/V* rv* I r 

d2 (3.24) 



f 2 

-{p-ki) [ki-{p-ki)] ■^P^pim,m)-^F\p_k^\{r]i,ri2)^ . 



Finally we get 
with 



{Qa{ki,r],)Q4k2,v*))cc = {27rfS'\ki + k2) Pa{h)^^ {ccclnh + acc) (3.25) 



ccc = ^ [-14 + xl (-7 + x^) (2 + x^)] , (3.26) 
and acc can be found in Appendix C. 

3.1.2 Off-diagonal contributions 

In this work, we consider the one-loop two vertices contributions. Besides the loop contributions from the same 
vertices, there are contributions involving two different vertices, as shown in Fig. 5. As described in the previous 
section, their explicit expressions can be read from these Feynman-type diagrammatic representations. 



ki 4 k k2 ki i k ki ki 4 ^ k2 



r], 9 • r], rj, 9 mi], rj, 

p-ki p-k] p-k] 

Figure 5: Diagrammatic representations of off-diagonal contributions. From left to right are contributions 
from (ab), (ac) and (be) combinations respectively. We do not show the (ba), (ca) and (cb) contributions 
explicitly. 

Contributions from (ab)-term are 

{Qa{ki,r]^)Qa{k2,'n^))^f,i 

'^'^^ J "^"^^mmJ (2^ ^P'^P- ^i)] (3.27) 

d^ '^'^ ' ^ d^ ^''^ ' ''^ ^ d^ ^"'^ ' ^ d^ ^p-*'^ I ^''^ ' ^ ^ 
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And 

{Qa{kl,V*)Q<7{k2,V*))ab2 

-co loo '^^^ ^ y (4^ ^ ■ " '''^^ (3-28) 

^^^^l^'^l''^*) ^^^^l^'^^''^^) ^^|p-fei|(m,^2) 

The (ba) contributions can be also read from Fig.5 which we do not write down exphcitly. Finally, we denote 
the contributions from (ab) and (ba) combinations as: 

{Qa{kl,r)^)Qaik2,V*))ab = i^Trf^Hkl + k2) Paih)^-^ {Cablnki + Q^fc) (3.29) 



And 



with 

_ -28x1 + 64 - 23x2 - 23 
~ 3840 

Contributions from (bc)-term are: 

{Qa{kl,r]*)Qa{k2,V*))bcl 



aril a'72 

And 

{QAkl,V*)Qa{k2,V*))bc2 

1 f cfip 



X -;f-Gki{vi,v*) GkA'n*,m) Fp{m,V2) ^F\p-ki\ijii,m)' 

arji ar]2 



(3.30) 



with „ . 

_ 44 - 104 - lla:^ + 49 

"""^ ~ 3840 ■ 

Contributions from (ac)-term are: 

{Qa{ki,v*)QAk2,r)*))aci 

dm / dm — / (p ■ fei) .3 

X :^GkAv*,m)GkAv*,V2) -^Fp{m,m) . ^ . ^|p- fei i , ??2 ) ■ 
am am aviam 

{Qa{kl,m)Qaik2,V*))ac2 

d d d^ 

Similarly, we denote the contributions from (ac) and (ca) as 

jj2 

{Qaiki,m)Qaik2,V*))ac = {^irf 5^ {ki + k2) Pc{ki)^^ (c„c In /ci + a„c) (3.33) 



(3.34) 



dm / dm — / 7^ [P • (P - fei)] (p • fei) (3 35) 



= -(27r)353(^^ + fe2) 85'cf £^ dm £^ '^'^''^J -^"^-^P- ^^^^ ■ (3.36) 
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While the total contributions from (be) and (cb) terms are, 

{Qa{ki,r]^)Qaik2,r]*))hc = i2nfd^{ki + fes) Pa{ki)^^ (c^c In fei + a^c) (3.37) 

with 



Che 



12x1 + 66a;f + 147^2 + 227 



(3.38) 



3840 
3,2 Final Results 

Collect (3.17)-(3.37) together, the final one-loop contributions from entropy mode to the adiabatic power spec- 
trum is ^ 

(Qa(?7*,fei)Qa(??*,fc2))i_ioop = (27r)=' 5=^(^1 + fe2)Pa(fei)^(cln/ci + a) , (3.39) 

TT €Cg 

with 

c = ^ [iSxt + bxl - 35) . (3.40) 

Note that the above result should be supplemented with a scheme-dependent UV renormalization constant. 
This result is consistent with previous analysis [50, 51, 52, 53, 54, 55, 56, 57, 58, 59]. In this work, the loop 
corrections show no IR divergences. This is mainly due to that in <C 1 limit, the leading-order interaction 
vertices are dominated by "derivatively-coupled" ones, which enhance the IR safety. Moreover, our result 
agrees with the bound obtained in [76, 77, 78], where perturbation theory breaks down. 



4. Conclusion and Discussion 

In this work, we calculate the one-loop corrections to the power spectrum of adiabatic mode from three-point 
adiabatic/entropy cross-interactions, at the time of horizon exiting. As mentioned before, in Cg ^ 1 models, 
there is enhancement of the effective couplings for interaction vertices, thus contributions (higher-order corre- 
lations functions and their loop-corrections) from multi-vertices diagrams are in general cannot be neglected. 
Precisely, one may expected the two-vertices one-loop contributions and the one-vertex one-loop contributions 
are of the same order. We find that, as the enhancement of non-gaussianities in models with non-canonical 
kinetic terms due to the enhancement of small speed of sound Cg, the loop corrections are also enhanced by 
both the slow-roll parameter e and Cg, more precisely, of order ~ Thus in general, the loop corrections 
can be large. Moreover, we find that the loop contributions are IR finite. Although in this work we shown 
these "large" loop effect through explicit calculations, they can be expected by a simple counting of the order 
of magnitude of the interaction vertices. The point is that, in relatively large parameter region we need to take 
into account the higher order loop corrections when we discuss models with small c^, which may be important 
in the investigation of large non-Gaussianities. 

There are several issues to be addressed. Firstly, the perturbations of inflaton fields and thus the correlation 
functions themselves are not observable directly. What is responsible for the CMB anisotropics and LSS is the 
curvature perturbation C,. Thus we should combine the perturbations of inflaton fields to yield the perturbation 
of curvature perturbation long after horizon exiting, e.g. by using 6N formalism [73, 74, 75], as in [55]. 
However, it is well-known that in multiple field inflation models, the curvature perturbation is not conserved on 
superhorizon scales, or in other words, there are also cross-correlations among different modes after horizon 
exiting. Due to these reasons, the complete analysis of loop corrections to the curvature perturbation in multi- 
field models is a complicate task which needs more subsequent work. 

In this work, we focus on the three-point adiabatic/entroy cross-interactions. While from (2.14), in general, 
there are adiabatic mode self-interactions. Moreover, in multiple-DBI models, these various interaction- vertices 
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give the same order of contributions, both to non-gaussianties and loop corrections. A complete analysis of 
loop corrections including contributions from adiabatic loops and also the relations with observables (e.g. the 
curvature perturbation and the CMB anisotropics) are also needed. 
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A. A Brief Review of "In-in" Formalism 
A.l Preliminaries 

The "in-in formalism" (also dubbed as "Schwinger-Keldysh formalism", or "Closed-time path formalism") 
[70, 71, 72] is a perturbative approach for solving the evolution of expectation values over a finite time interval. 
It is therefore ideally suited not only to backgrounds which do not admit an S-matrix description, such as 
inflationary backgrounds. 

In QFT apphed in the calculation of S-matrix in particle physics, the goal is to determine the amplitude for 
a state in the far past to become some state |V'') in the far future. 



Here, conditions are imposed on the fields at both very early and very late times. This can be done because that 
in Minkowski spacetime, the states are assumed to be non-interacting at far past and at far future, and thus are 
usually taken to be the free vacuum, i.e., the vacuum of the free Hamiltonian Hq. The free vacuum are assumed 
to be in "one-to-one" correspondence with the true vacuum of the whole interacting theory, as we adiabatically 
turn on and turn off the interactions between t = —oo and t = +oo. 

While the physical situation we are considering here is quite different. Instead of specifying the asymptotic 
conditions both in the far past and far future, we develop a given state /orwarJ in time from a specified initial 
time, which can be chosen as the beginning of inflation. In cosmological context, the initial state are usually 
chosen as free vacuum, such as Bunch-Davis vacuum, since at very early times when perturbation modes are 
deep inside the Hubble horizon, according to the equivalence principle, the interaction-picture fields should 
have the same form as in Minkowski spacetime. 

A.2 "In" vacuum 

The Hamiltonian can be split into a free part and an interacting part: H = Hq+Hi. The time-evolution operator 
in the interacting picture is well-known 



where subscript "I" denotes interaction-picture quantities, T is the time-ordering operator. Our present goal is 
to relate the interacting vacuum at arbitrary time \^i{t)) with free vacuum |0i) (e.g., Bunch-Davis vacuum). 



(V^'|5|V) = (V^'(+oo)|V(-oo)). 




(A.1) 
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The trick is standard. First we may expand \Qi{r])) in terms of eigenstates of free Hamiltonian Hq, |r2i(r/)) = 
\ni) {ni\Qi{r])), then we evolve \Qi{r])) by using (A.l) 

\nii7]2)) =U{v2,m)mm)) = |0i)(0i|f^i) + ^e+^^"(''2-''i)|ni)(ni|J^i(7yi)). (A.2) 

n>l 

From (A.2), we can immediately see that, if we choose r/2 = — oo(l — ie), all excited states in (A.2) are 
suppressed. Thus we relate interacting vacuum at r/ = — oo(l — ie) with free vacuum |0) as 



|J^l(-oo(l-ie))) = |0i) (Oi|17i) 
Thus, the interacting vacuum at arbitrary time r] is given by 

|VAC,in) = \ni{v)) = U{ri, -oo{l - ie))\Qi{-oo{l - ie))) 



(A.3) 



T exp 



oo(l— ie) 



drj'Hiiirj') |0i)(0i|f^i) 



(A.4) 



A.3 Expectation values in "in-in" formalism 

The expectation value of operator 0{7]) at arbitrary time rj is evaluated as 
(VAC, in 1 0(r/) I VAC, in) 



(VAC, in I VAC, in) 



Oi 



Texp i 



oo(l+je) 



dr,'Hii{r^') Oi(??)Texp 



oo(l— ie) 



dr]'H,i{r,') 



(A.5) 



Oi 



where T is the anti-time-ordering operator. 

For simplicity, we denote —00(1 — «e) = — cxd+ and — cxd(1 + ie) = —cxd^, since, e.g., —cxd+ has a positive 
imaginary part. Now let us focus on the time-order in (A.5). In standard S-matrix calculations, operators 
between (0| and |0) are automatically time-ordered. While in (A.5), from right to left, the time starts from 
infinite past, or — cxd+ precisely, to some arbitrary time rj when the expectation value is evaluated, then back 
to —00^ again. This time-contour, which is shown in Fig. 6, forms a closed-time path, so "in-in" formalism is 
sometimes called "closed-time path" (CTP) formalism. 



Im 7} 



-oo(l-ie) 



-oo(l+ie) 

Figure 6: Closed-time path in "in-in" formalism. 



Re?/ 



A.4 Perturbation theory 

The starting point of perturbation theory is the free theory two-point correlation functions. In canonical quan- 
tization procedure, we write a scalar field as 

(l)k{'n) =u{k,'q)ak + u*{k,r])ai^, (A.6) 

where u{k,rj) is the mode function for (pkiv) (i^i practice, Uk{i]) and u^(t/) are two linear-independent solu- 
tions of equation of motion for 4>k{r]), which are Wroskian normalized and satisfy some initial or asymptotic 
conditions ). 
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The free two-point function takes the form 

{0\(pk,{m)(PkM\0) = {2Trf5^{ki + k2)Gk,{vi,m), 

with 

In this work, we take (A.7) and (A. 8) as the starting point. 
Now Taylor expansion of (A.5) gives 

• Oth-order 



(A.7) 
(A.8) 



Ist-order (one interaction vertex) 

, (1) 

2nd-order (two interaction vetices) 

(2) 



2 Re 



(0) 



-oc 



(oiiai(77)|oi). 



dry' (Oi di{t)H,i{t'] 



Oi 



(A.9) 



(A. 10) 



-2 Re 



r dri' r dv" (Oi di{v)Hii{r)')Hii{ri") 

J-00+ J-00+ \ 

drj' r drj" (Oi H-a{r]')di{r])Hii{r]" 



Oi 
Oi 



(A. 11) 



Here in this work, since we are considering one-loop contributions from three-point interaction vertices, (A. 1 1) 
is needed. 

A.5 Positive/negative-fields notation 

There is a formulation of "in-in formalism" in terms of "doubled fields" in the literatures, mostly applied in the 
path-integral formulation. The reason is that, the closed-time path (CTP) is imagined to circle the time axis in 
the complex 77-plane (see fig.6). Accordingly, we distinguish operator values from the upper and lower branches 
by labeling a "+" on the upper increasing-time contour and a "— " on the lower decreasing-time contour. Then 
it is convenient to introduce a new "time-path-ordering operator" Tp, which orders operator along the "CTP- 
path", as shown in Fig.6. Obviously, Tp acts on " + " fields as time-ordering operator T, and acts "— " fields as 
anti-time-ordering operator T. Thus (A.5) can be recast into a more convenient and familiar form 



Tpexp(^i j' dv'Hr{r,')^ Oi(r?) Tp exp (^-i J'' d^'H+in'] 







Tp di{r]) exp 



dv' [H+irj')-Hr{r^')] 



(A. 12) 



Now, since we have two types of fields, the contractions between different pairs yields four kinds of correlation 
functions. 



T^pKi^i)KiV2) 0) ^ (27r)^5^(fci +fe2)G±±(77i,r72), 



where 



G'^~{m,V2) = G^{r]i,r]2)9{r]i - m) + Gl{r]i,7^2)d{m - m] 
G^~im,V2) = Gfc (r/i,r?2), 

G^^{m,m) = Gl{r]i,m), 



(A. 13) 



(A. 14) 
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and 

Gk{vi,V2) = Uk{vi)uUV2) , ^^^^^ 

Gtim^m) = u*kim)ukim) = {Glim^m))* = Gl{m,vi) , 

are the vacuum Wightman functions. Now the perturbation calculations are standard, after expanding the Tp- 
ordered operator exponential in (A. 12). 

The "doubled-field" notation is a convenient "book-marker" in path-integral formulation of "in-in formal- 
ism" in order to put everything into "one" exponential, while this notation itself is not necessary in principle. 
Especially, in operator formaUsm, there is no need to take this notation at all (see also [58] for a discussion on 
this point). 



B. Momentum Integral 



In this work, the momentum loop-integrals have the structure as 



^f{k,p, \p-k\) 



(B.l) 



where k is the momentum of external line, and p is the loop momentum. In this appendix, we briefly describe 

the methodology for evaluating such type of integrals. 

In general, one may use spherical coordinates {p, 0, (p) to calculate (B.l). However, it is more convenient 
to introduce a new variable 

q = \p — k\ 



2kp cos 6 , 



(B.2) 



and change coordinate {p, 6, 0) into {p, q, (j)). Note that this can be done only when / 0. The integral measure 
is simply d^p = ^dpdq d(j). Thus, the original integral (B.l) becomes 

+00 rp+k 



1 



The integral region for p and q is shown in Fig.7. 



dp I dqpqf{k,p,q) 

l\p-k\ 



(B.3) 




Figure 7: Integral region for variables p and q. 



One may use (B.3) as the starting point. However, it can be seen directly from Fig.7 that sometimes it is 
more convenient to use "rotated" variables 



x = q + p 
1 

P = o 



(x-y), 



y = q-p, 
1 

^ = n 



x + y), 



(B.4) 
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with integral measure dpdq = ^dxdy. Finally, the integral can be put into the form 

^^^^ = 8(2^ r '^'^ f-k ^""^ ~ ^ ^^"^ ~ ^^"^ ^ ^0 ■ ^^'^^ 

C. Finite parts 

For clarity, we collect the finite constant parts of the loop integrals. 

-59 - 256 In 2 - (59 + 256 In 2)xl + (464 - 1600 In 2)xt + (-30 + 64 In 2)xl 



an 



abb 



an 



(Xab 



a„ 



122880 

(l + x^) (331 + 30(-21 + 40 In 2)xt) 

460800 ■ 
439 - 3840 In 2 + (439 - 3840 In 2)xl - 120(-37 + 160 In 2)xt + 90(-33 + 64 In 2)xl 

1843200 ' 
-2687 - 43680 In 2 + (34913 + 4320 In 2)xl + 360(-47 + 40 In 2)xt + (1410 - 1280 In 2)x^ 

614400 

5647 - 6720 In 2 + ar^ (5647 - 6720 In 2 + 60^^ (23 - 40 In 2 + 8 In 2x1) ) 

230400 ' 
xl {30x1 (^2(11 + 32 In 2) -452 -80 In 2) - 11161 - 2640 In 2) + 4439 + 11760 In 2 

921600 ' 
-120x^(28 In 2 - 11) + 20x^(21 + 36 In 2) + x^(601 - 2760 In 2) + 601 - 2760 In 2 

460800 ' 
120x^ (12 In 2 - 19) + 60x^(77 + 132 In 2) + x^ (17640 In 2 - 4849) - 1249 + 27240 In 2 



460800 



(C.l) 



xl (-6x2 ^^2 + 22-18 In 2) + 347 + 30 In 2) + 137 - 210 In 2 



a = 

5760 



"be 

And 
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